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ABSTRACT 

Using the MIPS instrument on the Spitzer telescope, we have searched for 
infrared excesses around a sample of 82 stars, mostly F, G, and K main-sequence 
field stars, along with a small number of nearby M stars. These stars were selected 
for their suitability for future observations by a variety of planet-finding tech- 
niques. These observations provide information on the asteroidal and cometary 
material orbiting these stars - data that can be correlated with any planets that 
may eventually be found. We have found significant excess 70 fim emission to- 
ward 12 stars. Combined with an earlier study, we find an overall 70 /zm excess 
detection rate of 13 ± 3% for mature cool stars. Unlike the trend for planets to 
be found preferentially toward stars with high metallicity, the incidence of debris 
disks is uncorrelated with metallicity. By newly identifying 4 of these stars as 
having weak 24 /im excesses (fluxes ~10% above the stellar photosphere), we 
confirm a trend found in earlier studies wherein a weak 24 /im excess is associ- 
ated with a strong 70 /im excess. Interestingly, we find no evidence for debris 
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disks around 23 stars cooler than Kl, a result that is bolstered by a lack of excess 
around any of the 38 K1-M6 stars in 2 companion surveys. One motivation for 
this study is the fact that strong zodiacal emission can make it hard or impossi- 
ble to detect planets directly with future observatories like the Terrestrial Planet 
Finder (TPF). The observations reported here exclude a few stars with very high 
levels of emission, > 1,000 times the emission of our zodiacal cloud, from direct 
planet searches. For the remainder of the sample, we set relatively high limits on 
dust emission from asteroid belt counterparts. 

Subject headings: infrared: stars — circumstellar matter — Kuiper Belt 



Introduction 



A planetary system is characterized by the properties of its parent star, by the num- 
ber and nature of its gas-giant and rocky planets, by the extent of its Kuiper and asteroid 
belts, and by the populations of gas and dust orbiting the central star. In the coming 
decade, astronomers will use a variety of techniques to address all these aspects of neigh- 
boring solar systems. Initial results for gas-giant planets are based on ground-based radial 
velocity searches. Eventually, nearby stars will be the targets of indirect and ultimately 
direct searches for terrestrial planets with the Space Interferometer Mission (SIM) and the 
Terrestrial Planet Finder (TPF). The Spitzer telescope is uniquely positioned to characterize 
the evolution, amount, structure and composition of the dust associated with Kuiper and 
asteroid belts around many types of stars, including those with and without the gas giant 
planets now being detected by the radial velocity technique. Guaranteed Time Obs erver 



(GTO) studies such as the FGK sample (IBeichman et al. 



2005M iBrvden et al. 



the Nearby Stars pro gram iGautier et alil2006h . plus the FEPS Legacy project (IMever et al.l 
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2004 ; iKim et al.ll2005l ). have conducted photometric surveys of about 200 nearby stars at 24 
and 70 /mi. The photometric survey discussed here uses Spitzer images at 24 and 70 /im 
to look for debris disks around an additional 82 stars, rounding out existing surveys of the 
closest stars. 

As the Spitzer programs are completed, we will be able to carry out statistical in- 
vestigations of the debris disk phenomenon in terms of the age, metallicity, and spectral 
type of parent stars. In particular, by nearly doubling the size of the existing sample of 
stars (relative to the on-going GTO/Legacy programs) we can hope to identify and im- 
prove the statistics of types of excess that appear to be rare based on existing IRAS or 
ISO observations, e.g . hot dust, extremely large disk to star luminosity ratios (Lri„ st ,/L* ) 
around mature stars (IFajardo-Acosta et al.ll2000l ; lHabing et al.ll200ll ; ISpangler et al.ll200ll ). 
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or low mass stars. The incidence of excesses at t he IRAS/ISO sensitivity level is about 15% 
( IBackman &: Parescd Il993l ; iBryden et al.l l2006af ) so that a total Spitzer sample of 250-300 
stars can hope to identify over 50 stars with excesses suitable for fut ure study. Stars with 



hot excesses (peak wavelength < 24 /xm) are con siderably rarer, 2-3% (IFajardo-Acosta et al. 



2000l ; lLaureijs et al.ll2002l ; iBeichman et al.ll2006l ). so that a survey of a large number of stars 
is needed to generate a statistically meaningful sample. 

As the statistics of planets build up, we will be able to correlate the properties of 
debris disks (total mass, physical configuration, composition) with properties of number, 
location, and mass of planets. Much lower dust masses can be de tected with Spitzer than 
was previously pos s ible, p articularly for solar-type and cooler stars. IBeichman et al.l (l2005bl ) 
and IBryden et al.l (j2006al ) have shown that with Spitzer instruments, we can reach just a 
few times the fractional l uminosity p redicted for our own Kuiper Belt (0.3 — 5 x 10~ 6 ; 
Backman fc Parescd 1 19931 ; ISternlll996l ). Determining how many mature stars like the Sun 
have Kuiper Belts comparable to our own is an impor tant ingredient in understan ding the 
formation and evolution of solar systems like our own (ILevison fc Morbidellill2003l ). 



Finally the Spitzer data will help us to understand the potential influence of zodiacal 
emission on the eventual direct detectability of planets. As hi ghlighted in a numb er of TPF 
studies, including the Precursor Science Roadmap for TPF (ILawson et al.l 12004 ) . the level 
of exo-zodiacal emission can affect the ability of TPF to detect planets directly, particularly 
for extreme cases with much greater dust contamination than in the Solar System. A com- 
plete census of potential TPF stars will assist in the eventual selection of TPF targets by 
determining or setting a limit to the amount of exo-zodiacal emission around each star. 

This paper focuses on the results of the 24 /xm and 7 /im survey using the Multiband 
Imaging Photometer (MIPS) instrument on Spitzer (MIPS; lRieke et al.ll2004l ). After describ- 
ing our target selection (§2]), we present these MIPS observations in §3J As will be discussed 
below, a number of sources in this sample appear t o be extended . These are highlighted in 
§31 but are discussed in detail in a separate paper (IBryden et al.l l2006bl ) . Follow-up obser- 
vations of sources with excesses using the IRS spectrometer are just now being completed; 
these will also be detailed in a later paper. In §Hwe discuss how our MIPS flux measurements 
constrain the dust properties in ea ch system. Fo r the sy stems identified here as having IR 
excess, combined with those from IBryden et al.l (l2006al ). §5] attempts to find correlations 
between the dust emission and system parameters such as stellar metallicity, spectral type, 
and age. Finally, in £0 we assess the influence of debris disks on the detectability of planets. 



-4- 



2. Stellar Sample 



Our sample is based on work carried out by radial velocity search teams (e.g. iMarcy et al. 



20041 ) and by the SIM and TPF Science Teams to identify the most suitable targets for the 
indirect or direct detection of terrestrial-mass planets (1-10 M e ). One target list consists 
of the 100 stars in the SIM Tier-1 sample which will be the most intensively o bserved stars 



in the two SIM projec ts dedicated to finding planets around nearby stars (IMarcy et al. 



2002; Shao et al. 



iooih . Q s ince the absolute astrometric signal from a planet scales as 
3 /zsec (At/pc) -1 x (api ane t/AU) x (Mpi an et/M e ) x (M /M*), the SIM teams are concentrating 
on the some of the closest, lower mass stars for their deepest surveys for terrestrial planets. 
Thus, the SIM list includes a number of late K and M stars not included in the other Spitzer 
samples or in TPF lists. 

We also draw from a number of lists prepared by Science Working Groups for the TPF- 
Coronagraph (TPF-C) and TPF-Interferometer (TPF-I) missions. Although the TPF lists 
are n ot definitive given the indeterminate status of the project, the outline of the sample is 
clear (jBeichman et al.ll2005d ; iTraub et al. 1120061 ). We start with F0 - M5 stars of luminosity 
classes IV or V and refine the list by making a few simple assumptions about the nature of 
planetary systems and the properties of TPF. Specifically, we 1) exclude stars with binary 
companions within 100 AU as being inimical to the formation or sta ble evolution of plan etary 
systems; 2) require that the angular extent of the habitable zone feasting et al.lll993l ); ~1 
AU for a 1 L Q luminosity star, scaled by the square root of the stellar luminosity) exceed 50 
milliarcsec; and 3) impose an outer distance cutoff of 25 pc (although we allowed a few F0-F5 
stars at distances as great as 30 pc to bring up their numbers). To enable good measurements 
with Spitzer we rejected stars with high levels of stellar and/or cirrus confusion based on 
examination of IRAS maps. 

Comparison of potential SIM and TPF targets in cirrus-free sky with the Spitzer Re- 
served Object Catalog (as of November 2003) showed 81 stars with spectral types ranging 
from F0 to M3.5, as listed in Tabled! One more star, GL 436, was added through a Director's 
Discretionary Tim e proposal after the discovery of a planet in this system was announced 
( IButler et al.ll20041 ). Binary companions within the 82 fields of view have also been included 
as secondary targets; six such companions are identified as bright enough for clear detection 
in both the 24 and 70 /im images^ The divergent proper motion of HD 48682B, an M0 star 



1 The merged, high priority target list for these projects is available at 
http: / / astrori.berkeley.edu/^gmarcy/sim_draft.html| 

2 The 24 /im image of HD 265866 has what appears to be an equal-brightness binary companion located 
40" NW of the target primary star. However, there is no visible or near-IR neighboring source. In fact, 
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30" to the NE of HD 48682, rules out a physical association between the two stars. Thus, 
HD 48682B is not included in this sample. The true binarity of the other six neighboring 
sources is verified via their Hipparchos distances and space motion measurements. Angular 
separations in these systems range from 10" to 100", with projected orbital separations be- 
tween 100 and 1000 AU. Data for the companions are listed separately at the end of Tables [1] 
and [21 With their inclusion, our total sample contains 88 stars within 82 targeted fields. 

Binned by spectral type, the SIM/TPF sample consists of 37 F stars, 19 G stars, 24 K 
stars, and 8 M stars. Typical distances range from 10 to 20 pc, closer for M and K stars and 
farther for earlier spectral types. Figure [1] shows the overall distribution of observed spectral 
types. Some basic parameters of the sample stars are listed in Table [H most importantly 
age and metallicity, which are also shown as histograms in Figures [2] and [21 There is no 
explicit target selection based on stellar age or metallicity, but known planet-bearing stars 
have been specifically included in a couple cases. Of this sample, only two stars (GJ 436 
and HD 147513) are already known to have planets; most of the other stars with planets are 
either too fain t, lie in cirrus-contamin ated regions, or are already observed in other Spitzer 



programs (e.g. 



Beichman et al.ll2005bl ). 



In this paper we first discuss t he 88 primary and secondary stars and then add in 
the stars observed in Bryden et al. (|2006a ) to increase the size of the sample for some of 



statistical discussions. 



3. Spitzer Observations 

All stars were observed with MIPS at 24 //m and, with one exception (HD 265866), at 
70 /im. In order to help pin down their stellar photospheres, four M stars - GJ 908, HD 36395, 
HD 191849, and HD 265866 - were also observed with the IRAC camera in subarray mode 
at all four of its wavelengths (3.6, 4.5, 5.8, and 8.0 /xm). Seven stars identified as having 
IR excess were observed with IRS, the Spitzer spectrograph, as follow-up observations, as 
detailed in a future paper. 



the second 24 /im source is a chance alignment with a passing asteroid. Software specifically dev eloped for 
locating asteroids relative to the Spitzer observatory (part of the Horizons package; lGiorginill2005f) iden tifies 
this object as asteroid #11847 ( "Winckelmann" : H=13.4, a=2.67 AU, e=0.065, i=10.23: lBowel]||l996h . 



- 6- 



3.1. Data Reduction 



3.1.1. MIPS Observations 



Overall, our data an a lysis is similar to that previously described in iBeichman et al. 



(2005b 


) and 


Brvden et al. 


(2006a) 



the DAT software developed by the MIPS instrument team (IGordon et al.ll2005l ). At 70 /xm, 
images were processed beyond the standard DAT software to correct for time-dependent 
transients, correction s which can significantly improve the sensitivity of the measurements 



( IGordon et al.l 12004 ) . For both wavelengths, aperture photometry was performed using 



apertures sizes, backgro und annuli, aperture corrections, and instrument calibration as in 



Beichman et al.l (j2005bl ). We find that the target lo cations in the 24 u m images are con 



sistent with the telescope pointing accuracy of <1" (IWerner et al.l 120041 ) . As such, we use 
the 24 /im centroid as the target coordinates for both wavelengths. Special consideration is 
made for the 6 resolved binaries in our sample. Instead of our standard method of aperture 
photometry with a surrounding sky annulus, the emission at the two stars' locations is fit 
with the instrument's point spread function (PSF). We find that for binary stars with small 
angular separations, simultaneously fitting of their overlapping PSFs results in much im- 
proved photometric accuracy. The agreement between PSF fitting and aperture photometr y 



(with appropriate aperture correction) for isolated stars is excellent (IGordon et al.ll2005l ). 



For all of the stars, the MIPS flux and noise measurements are listed in Table |2j 



3.1.2. IRAC Observations 



The IRAC sub-array images of t he four M stars were r educed following the technique 
described by the FEPS Legacy team ([Carpenter et al.l 120061 ) . The pixel sizes are corrected 
for distortion and a pixel-phase correction is made to channel 1. Stellar fluxes are measured 
within an aperture of 10 pixels (=12"), with a background annulus from 10 to 20 pixels. The 
photometric measurements each star at the four IRAC wavelengths are listed in Table 



3.2. Photospheric Extrapolations and Limits on 24 |im Excess 

To determine whether any of our target stars have an IR excess, we compare the mea- 
sured photometry against predicted photospheric levels. A detailed descriptio n of our stellar 



atmosp here fitting, as applied to F5-K5 stars, is presented in the appendix of iBryden et al. 



( ]2006al ). The stars observed here, however, span a greater range of spectral types than 
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previously considered. In particular, our sample contains late K and M t ype stars with 
numerous broad molecular features for which the stellar models (IKuruczl 120031 ) begin to lose 
their accuracy. 

This accuracy can be directly assessed by examining how well the observed flux levels 
match those predicted. We use the ratio, Fmtps/-P *-24ahii, to as s ess the photospheric ex- 
trapol ation using the fact previously established in iBryden et al.l (j2006al ). iBeichman et al. 
( 120061 ). and earlier references cited therein, that excesses at 24 /im are rare (~1%). Fig- 
ure |H shows the distribution of this ratio. After excluding one outlying star with a strong 
24 /im and 70 /im excess (HD 109085), the 88 flux measurements at 24 /im have an average 
Fyiips I of 1.01. The dispersion of Fmips/ in Figure H] is 0.10, which is rela - 
tively large compared to the previous result for just F5-K5 stars (0.07; IBryden et al.ll2006al ). 
We identify three causes for this larger dispersion: 



3.2.1. Quality of Near- IR Photometry 

The SIM/TPF sample contains a number of nearby stars that are brighter than the 
stars in the FGK survey. Stars brighter than about K s = 4 mag have saturated 2MASS 
measurements resulting in large photometric uncertainties (~0.25 mag). For several of these 
stars, particularly the early F type stars which are the brightest in the sample, Johnson 
.fT-band photometry is available in the literature with much better accuracy (~0.05 mag) 
than the saturated 2MASS values, but worse than the best 2MASS values (~0.03 mag). 
In such cases, the saturated 2MASS values are supplanted by the better data. However, 
the uncertainty in the near-IR photometry for the remaining 2MASS-saturated stars causes 
difficulty in extrapolating to longer wavelengths and results in a greater dispersion than 
when only stars with high quality 2MASS data are used, ct(Fmips/-F*-24/ziii) = 0.09 vs. 0.07 
for types F5-K5. 



3.2.2. Intrinsic Variability 

Stellar variability between the epochs of the Spitzer data and the photometry used to 
estimate the photospheric contribution could account for some of the dispersion in FmipsI F*~ 
24/xm. To investigate this possibility we examined the Hipparcos photometry for the 66 stars 
of our sample for which these data are available (ESA 1997). Only three stars (GL 436, HD 
79211 and HD 265866) showed a scatter in mag(Hp) in excess of 0.02 mag while the vast 
majority had scatter less than 0.01 mag. In particular, none of the ten stars younger than 1 
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Gyr and thus possibly more variable than the rest of the sample, showed variability above 
this level. GL 436 and HD 265866 are faint V~10 mag M stars so that the level of Hipparcos 
scatter is not significant. The large scatter for HD 79211 (0.23 mag) is due to multiplicity 
and is also not significant. As discussed below, one non-Hipparcos star, HD 38392, shows a 
low level of variability and a correspondingly larger deviation in Fuips/ -F*-24/im. 



3.2.3. Quality of Photo spheric Models 



The ability to extrapolate from visible and near-IR photometry to MIPS wavelengths 
appears to be an issue for spectral types later than the F5-K5 range used in the FGK survey. 
For all stars with accurate 2MASS data, Figure [5] plots the directly observable K s — [24] color, 
a quantity independent of the stellar atmosphere models. With the exception of the M stars, 
all of the averages are consistent with a constant color of K s — [24] ~ 0.02 ± 0.02. Most 
interestingly, an apparently abrupt transition occurs between the late K stars and M stars, 
with the average K s — [24] color jumping up ~0.4 magnitudes for the cooler stars. T his trend 
of redder K s — [24] for later spectral types was first noticed by lGautier et al.l (120061 ) whose M 
star data are shown for comparison. 

We next consider the ratio of the observed flux at 24 fim to that predicted by pho- 
tospheric models (Fmips/ -F*-24/im) as a function of spectral type. Solar-like stars (types 
F5-K4) have an overall average of FMips/-F*-24/im = 0.98 ± 0.01 with a dispersion of 5% 
among the stars with good 2MASS data and excluding stars with excess emission at 70 /zm 
(identified in §3.3p . For F0-F4 stars, the observed fluxes are marginally higher than those 
predicted, with an average i*MiPs/-^* _ 24//m = 1.03 ± 0.02. For late K stars with good 
2MASS observations the observed fluxes are consistently below expectation with an average 
F MIPS /F*-24/mi of 0.87 ± 0.03. Since the observed K-[24] color is flat (Fig. EJ), this offset 
is likely a fault of the photospheric modeling or of our fitting procedure. Not surprisingly, 
the models have the greatest difficulty with the M stars which have average i 7 Mips/-F 1 *-24/zm 
= 1.16 ± 0.06. This dif ference between measu red and predicted fluxes for the M stars re- 
mains even if Next Gen (IHauschildt et al.lll999l ) models are used instead of Kurucz models. 
However, with an accurate determinat ion of each star's effective t emperature and with more 
advanced stellar models fPHOENIX: iBrott fc Hauschildtl I2OO5I ) . bautier et all J2006h were 
able to fit the 24 /zm colors of M stars. 

Thus, knowing that these trends in i 7 Mips/-F 1 *-24/im exist and may ultimately be ex- 
plained with better modeling, we can compare each star with the average K s -[24] color 
within its spectral type bin (Fig. [5]) to look for dust excesses. With this methodology, we 
find no evidence for a 24 /im excess toward any of our M stars. This negative result is 
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strengthened when IRAC photometry is available. For the 4 M stars with IRAC data (Tar 
ble|3]), the inclusion of 3.5-8.0 /im fluxes into the fit modifies the average of -Fmips /F*-24/xm 
from 1.32 with a dispersion of 0.21 to i 7 Mips/-^*-24/im=0.92 with a dispersion 0.04, confirm- 
ing that the M stars' 24 /im fluxes are consistent with emission from the stellar photosphere 
alone. 



3.2.4- Presence of a Weak Excess at 24 fim 

Finally, the third reason for increased dispersion in the FMips/-F*-24/mi values, in ad- 
dition to poor near-IR photometry and less accurate stellar photospheres or model fitting 
for late type stars, is the presence of weak but real excess emission from dust toward some 
stars. In §3.31 we will identify some of our target stars as having strong excess emission at 
70 /an. Only one of these objects, HD 10908qj (the labeled value in Figure HJ) also has an 
immediately obvious IR excess at 24 //m. 

Taken in composite, however, the stars with 70 /im excesses tend to have a weak 24 /im 
excess. Considering only F0-K5 stars with good near-IR photometry, those with 70 /im 
excess have an average FMips/-^*-24/im value 9% higher than those without. A similar 
general trend was previously noticed by iBryden et al.l (j2006al ) and was confirmed in the 
IRS spectra of individual object s with 70 /im excess, which tend to rise above the stellar 
photosphere longward of 25 /tm (IBeichman et al.ll2006l ). Combining the F0-K5 stars in this 
sample with those from IBryden et al.l (l2006al ). we are no longer limited by small number 
statistics and the correlation between 70 /tm and 24 /tm excess becomes significant at the 
3-cr level. The average 24 /im excess for stars with 70 /tm excess is 0.079±0.026 times the 
stellar flux. 



To assess the significance of a possible 24 /im excess on a star-by-star basis, we define 
the parameter x 24 which corresponds to the n — a significance of any deviation from the 
expected photospheric value: 

F24 — F* 



X24 = 



1 



0~24 



where F 24 is the measured flux, is the expected stellar flux, and <t 24 is the noise level, all 



An excess was first detected around HD 109085 (= 77 Crv) by IRAS (IAumannl ll988; Stcnc el fc Backman 



19911 ) and subsequently with SCUBA at sub-mm wavelengths (jSheret et al 



200 



Wvatt et al 



2005F ". Con- 



sistency between the IRAS flux at 25 /im and the MIPS 24 /im flux measured here depends strongly on the 
application of color corrections which are functions of the assumed dust temperature. For dust temperatures 
around 200-400 K and assuming the photospheric value given in Table [21 we find good consistency between 
the two measurements. 
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at 24 /im. A similar definition follows for 70 jum (equation [2]). We take the noise to be the 
larger of either 4% for sources with good 2MASS or Johnson data or 8% for sources with 
poor near-IR photometry. These values are based on the dispersions in F M ips/-F 1 *-24/im for 
the stars without 70 /im excesses. Ignoring the previously discussed late K and M stars with 
poor photospheric extrapolations, we find that the deviations from photospheric emission 
skew sharply to positive values for stars with 70 /im excess (black shading in Fig. HJ). Using 
this analysis we identify statistically significant 24 /im excesses accompanying a stronger 
70 jum excess around two stars: HD 25998 (3.6 a) and HD 40136 (3.2 a) in addition to 
HD 109085 discussed earlier. At slightly lower significance we find hints of a 24 /im excess 
for HD 199260 (2.7 a) and HD 219482 (1.9 a) which the accompanying 70 /im excess suggests 
could be real. 

A number of other stars show strong deviations from photospheric values without an 
accompanying 70 /im excess: the deviant -Fmips /.F*-24/im values of the M stars have already 
been discussed and attributed to poor photospheric extrapolation; HD 38392 has an ap- 
parent 30% excess at 24 /im which we attribute to the difficulty of obtaining an accurate 
measurements due to a) saturated 2MASS measurements, b) proximity to a nearby, bright 
companion (HP 38393) and c) th e possible variability of the star itself at the 5% peak-to- 



peak level flNitschelm et all 120001 ) . Finally, HD 23249 and HD 55892 have 2 < ^24 < 3 and 
X7o < 2. These deviations could simply be statistical fluctuations or they could be hints of 
an excess like that seen toward HP 69830 which is prominent only in the 8-34 /im region 



but not at 70 /im (jBeichman et al.ll2005a ). Without additional data, e.g. IRS spectra, we 



cannot assess the reality of the excesses around these last two stars. 



3.3. Detection of 70 /im Excess 

Having used 24 /im fluxes to test the accuracy of our stellar photosphere predictions, 
we next consider the frequency and strength of excess emission at 70 /im. The distribution 
of 70 /im flux densities relative to the expected photospheric values is shown in Figure [6J 
Unlike the tight distribution of flux ratios at 24 /im, many stars have 70 /im flux densities 
much higher than expected from the stellar photosphere alone. In several cases, the flux is 
more than an order of magnitude greater than expectation. Twelve of these stars will be 
identified in the following as having statistically significant IR excess. Excluding these stars 
with excesses and those with signal-to-noise ratio (S/N) < 3, the average ratio of MIPS 
flux to predicted photosphere is |Fmips/-F*-70/zhi| = 1.02 ± 0.05, consistent with the overall 
calibration. 

The dispersion in the 70 /im data is ~40% (excluding the stars with excesses), consid- 
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erably higher than that in the 24 /im data. An analysis of the noise levels in each individual 
field is required to assess whether the IR excesses are statistically significant. Many contri- 
butions to the overall error budget must be considered including those arising from stellar 
photosphere modeling, instrument calibration, sky background variation, and photon de- 
tector noise. At 70 fim, the calibration uncertainty and the background noise within each 
image are considerably larger than at 24 /im. On top of an assumed calibration uncertainty 
of 15%, we directly measure the standard deviation of the background flux when each field 
is convolved with our chosen aperture size. This background noise, which ranges from ~2 
to 20 mJy with a median of 3.7 mJy, is due primarily to extragalactic source confusion 
and cirrus contamination, rather than photon noise, and hence cannot be greatly reduced 
by additional integr ation time (for a more detailed analysis of the 70 /im noise levels, see 



Bryden et al.ll2006al ). Based on this measured background noise, we determine the S/N for 
each star, as listed in Table [2J Despite the high level of noise in some fields due to cirrus 
contamination and/or background galaxies, 72 out of the 87 stars in our sample with 70 /im 
data are detected with signal-to-noise ratio greater than 3. The median S/N for all of our 
target stars is 6.6, excluding the sources identified as having excess emission (which have a 
median S/N of over 20). 

Adding both background noise and calibration error together gives us a total noise 
estimate for each 70 /im target. In Table [2] we list these noise levels, along with the measured 
and the photospheric fluxes, for each observed star. We use these noise estimates to calculate 
Xro which corresponds to the n — a significance of any deviation from the expected level of 
photospheric emission: 

F 70 - F+ 

X70 = (2) 

070 

where F 70 is the measured flux, is the expected stellar flux, and a 70 is the noise level, 
all at 70 /im. Based on this criterion, we find that 12 out of 88 stars have a 3-a or greater 
excess at 70 /im: HD 25998, HD 38858, HD 40136, HD 48682, HD 90089, HD 105211, 
HD 109085, HD 139664, HD 158633, HD 199260, HD 219482, and HD 219623. In a sample 
of 88 stars there should be fewer than 1 star with a spurious excess on purely statistical 
grounds. Although cirrus or extragalactic confusion could produce spurious excesses, careful 
examination of each of the 70 /im images suggests that 

this is unlikely in the vast majority of cases. For example, the 70 /zm emission is well 
centered on the 24 /im positions, typically within~l". A number of weak excesses could 
have escaped detection under these criteria. Observations at higher sensitivity or at higher 
spatial resolution will be needed to identify these. The detection rate o f 70 /im excess within 



this sample is 14 ± 4%; combined with the sample of iBryden et al.l (j2006al ). this gives an 
overall detection rate of 13 ± 3% for cool stars. 
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The eight newly discovered IR-excess stars mostly have 70 jttm fluxes less than 100 mJy, too 
dim to have been be detected by IRAS. The notable exception is HD 105211 which has a 
very strong 70 /im flux (~500 mJy), but lies near a bright infrared source (CL Cru); with 
a separation of 2.4', this source is easily resolved in the MIPS image (Figure Ej), but still 
contaminates the large IRAS beam. For the stars without any significant excess emission, 
3-cr upper limits on possible excess flux typically range from 0.2 to 1.0 times the stellar flux, 
with a median of upper limit of 0.6 F*. 

Although the telescope resolution at 70 /im is relatively poor (FWHM of ~17"), several 
of these sources appear to be slightly extended in t he MIPS images (ma rked with superscript 
/in Table [2]). As discussed in a separate paper (IBryden et al.ll2006bl ). examination of the 
images of these marginally resolved sources do not indicate contamination by background 
objects, e.g. cirrus or galaxies, but rather that the objects possess to be truly extended 
disks. For one of the five Spitzer-resolved disks, HD 139664, a Hubble Spac e Telescope 



imag e shows the same orientation of the disk in the visible as in the infrared (IKalas et al. 



20061 ). At distances of 10-20 pc, the resolved disks have apparent radii of 100's of AU. As 
discussed below, maintaining a temperature of ~50 K at these distances (warm enough to 
emit strongly at 70 /im), requires relatively small grains with low emissivities. 



4. Properties of the Detected Dust 



Beyond our initial goal of detecting IR excesses, we are interested in determining the 
properties of the dust in each system - its temperature, luminosity, mass, size distribution, 
composition, orbital location, etc. For the 12 stars with significant IR excess, Table H] lists 
the excess 70 /im emission and measurements of or limits to the 24 jum excess. Seven of 
the 12 excess sources have an excess measured only at one wavelength (70 /im), restricting 
our ability to place limits on key quantities. The observed flux can be translated into the 
total dust disk luminosity relative to its parent star only when some as sumption is made for 
the dust temperature (e.g. iBeichman et al.ll2005bl : IBryden et al.ll2006al ). The minimum disk 
luminosity as a function of 70 /im dust flux density is obtained by setting the emission peak 
at 70 /im (or, equivalently, setting Td us t = 52.5 K): 



L 



dust 



minimum) = 10 



_ 5 /5600 K 



Fj,(70 /xm, dust) 
F„(70 /xm,*) 



(3) 
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Based on this equation, a minimum L^t/L* is calculated for each of our target stars iden- 
tified as having an IR excess at only 70 /an (Table HJ). For the other IR-excess stars, we 
also have at least a rough (2-cr) detection of the 24 /zm excess, as discussed in the previous 
section. With two wavelengths of excess measurement, the dust emission can be fit with a 
representative dust temperature, otherwise only upper limits can be obtained (Table HI also 
see Figure [TT1 below) . Figure [8] shows a spectral energy distribution (SED) for two of these 
stars, HD 219482 and HD 40316, fit with temperatures of 170 and 80 K respectively. In the 
cases with a fit dust temperature, Table @] lists the ratio of the integral under the stellar (T*) 
and dust (Tdust) blackbodies as a proper, rather than minimum, estimate of Ldust/L*: 

'T dust \ 4 - 1 \ F u (70 fim, dust) 



^dust 



U \ % J \ e**-l J F u (70fim,*) 
where x = hvjkT = 205.7 K/T at 70 fim. 



(4) 



For each of the stars with excess emission (plus those from iBryden et aJj (j2006af ) , Fig- 
ure [9] shows the total dust area and radial location of the emitting material. Despite 
the expectation that only large grains should be seen around mature systems due to loss 
mechanisms such as Poynting-Robertson drag and radiation pressure, small grains must 
be considered as a ser ious possibility given their presence in a mature star like HD 69830 



(iBeichman et al.ll2006l ) and because of the large extent of at least some the dust disks. Thus, 
the orbital location of the emitting material can only be calculated if some assumption is 
made for the dust's emissivity; small grains are less efficient at emitting infrared radiation, 
resulting in a higher temperature for a given orbital radius. For a given dust temperature, 
the orbital radius decreases with emissivity as e~ 5 . As such, for each calculated dust tem- 
perature or upper limit in Table HI two locations are plotted - the location if the emitting 
material is large blackbodies (lower axis) and the location if it is small grains with emissivity 
= 0.01 (upper axis). 

The dust area and mass have a similar ambiguity based on the unknown dust size/emission 
properties. The dust area in Figure [9] (left axis) is calculated under the assumption of black- 
body grains (unity emissivity); a lower emissivity would increase the dust area in direct 
proportion to e. Lastly, the dust mass (right axis) is based on an assumed typical grain size 
of 10 fim. An unconstrained amount of mass is contained in the larger parent bodies whose 
collisions produce the emitting dust. 

For stars with no detected emission, 3-cr upper bounds on the 70 fim fluxes lead to 
upper limits on L dust /L+ as low as a few times 10~ 6 , assuming a dust temperature of ~50 K 
(Table [2]). Although we cannot rule out cold dust at >, 100 AU, we are placing constraints 
on dust at Kuiper Belt distances to ~10-100 times the level of dust in our solar system. 
The constraint on asteroid belt-type dust is less stringent, ~100-1000 times our zodiacal 



-14- 



emission. 



4.1. Comparison with Sub-mm Observations 



The dust properties can be further constrained with sub-mm flux measurements. When 
available, longer wavelength data can help place lower limits on dust temperature, upper 
limits on the dust luminosity, and, with some assumption on the grain emissivity, outer 
limits on the disk extent. For most of our stars with IR excess, large amounts of cold dust 
emitting at longer wavelengths cannot be ruled out, but three of the 70 /zm excess stars have 
been observed at 450 and 850 /iin with JCMT /SCUBA, with two detections (HD 48682 and 
HD 109085) and one upper limit (HD 139 664; Sheret et al.l 2004 ) . Combining their sub-mm 



data with the infrared fluxes from IRAS, Sheret et al. 



(12004h modeled the SEDs for these 



stars, obtaining dust temperatures of 99 and 85 K for HD 48682 and HD 109085, respectively. 



4.2. Comparison with Visible Observations 



A nearly edge-on disk around HD 139664 h as recently been d etected at visible wave- 
lengths using the Hubble Space Telescope (HST) (IKalas et al.ll2006l ). The HST image of the 
disk shows a well defined inner edge around 83 AU and an outer edge that extends out to 
about 109 AU. If the dust associated with the 70 /mi emission is located in this ring, then a 
very large surface area of 10 /im grains would be required for the IR emission, since dust at 
this distance would have an equilibrium temperatu re of 30-35 K using a standard radial power 
law for grain temperature (jBeichman et al.l 120061 ) . Smaller, 0.25 /im grains, on the other 
hand, give temperatures of 69 and 77 K at the ring boundaries (cf. our 3-sigma upper limit 
of 78 K; Table HJ). Using a simple relationship f or dust mass, Mduxt 



4 D2Fy(dust) 
^P a oram QabsD (Tdust i 



and standard silicate absorption efficiencies, Q a bs (IDraine fc Ledll984j ; lBeichman et al.ll2006l ). 
yields a mass of 2.4 x 10~ 3 M e in large grains or 1.6 x 10 _4 M e in small grains, where we 
have taken a grain density of p=3.3 gm cm~ 3 and a distance, D = 17.5 pc, for this star. 

The radiative blowout si z e for grains a r ound an F5 star is approximately 1 /im (§ V.B.I 



in 



Backman &: Parescd (119931 ): iBurns et al. 



of small grains seen toward Vega (ISu et al. 



1979k In contrast to the spherical distribution 
2005) which is probably due to a recent catas- 



trophic event, small grains should be quickly ejected from the presumably stable ring system 
of HD 139664. IRS spectroscopy and/or millimeter spectroscopy would help distinguish 
between the large and small grain models. 
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5. Correlation of Excess with System Parameters 



To understand the origin and evolution of infrared excess, we now consider the properties 
of the sample stars and how they correlate with excess detection. Specifically, we examine 
the correlation with three variables: metallicity, age, and spectral type. These parameters 
are list ed for each star in T able [TJ Where appropriate we merge the present sample with 
that of iBryden et al.l (|2006al ) to improve the significance of any statistical conclusions. 



5.1. Metallicity 



Table [T] lists the metallicity information obtained from the literature for each of our 
target stars (number of independent [Fe/H] estimates, their average, and their r.m.s. scatter). 
Figure [3] shows a histogram of these metallicity values, the majority of which range between 
- 0.4 to +0.1 dex with a mean value of somewhat below solar. The stars with IR excess are 
identified with vertical arrows with the length of each arrow proportional to the strength 
of 70 fim excess. We find no correlation between metallicity and IR excess. The average 
[Fe/H] is -0.15±0.03 for all the observed stars and -0.17±0.04 for the stars with excess - an 
insignificant difference. 



These data are combined with IBryden et al.l (j2006al ) sample to show the fractional 
incidence of disks as a function of [Fe/H] (Figure [T2|) . A x 2 test shows that the distribution 
of disks in our three metallicity bins with significant number of stars (-0.75 <[Fe/H]< 0.0) 
is indistinguishable from flat. The lack of correlation between IR excess and metallicity is 
in sharp contrast wit h the well known correlation betwee n extrasolar gas giant planets and 
host s tar metallicity (lGonzaleall997l ; ISantos et al.l 120011 ) . In particular, iFischer fc Valenti 
(120051 ) find that the probability of harboring a radial-velocity detected planet increases as 
the square of the metallicity (Figure IT21. Although one might suspect giant planets and 
debris disks to be related, a similarly strong correlation between dust and metallicity can be 
confidently ruled out. A \ 2 comparison between the disk and planet distributions in same 
3 [Fe/H] bins suggest that there is only a 0.3% probability of these being drawn from the 
same distribution. 

The lack of correlation is furth er confirmed via Mon te Carlo simulations. Again using 
our dataset combined with that of IBryden et al.l (l2006ah (giving a total of 19 excess stars 
out of 151 observed), the correlation coefficient, r, is calculated for 10,000 random samples 
of stars. The histogram of the resultant r values is shown in Figure H31 under two different 
assumptions. In one case, the stars with excess are chosen randomly (left histogram centered 
on r=0); in the other, the stars are chosen with weighting proportional to their metallicity 
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squared (right histogram with average r=0.33). The correlation coefficient observed within 
our data (vertical arrow) is inconsistent with the strong metallicity dependence observed 
within planet-bearing stars. 

This lack of correlation may reflect the different formation histories of giant planets and 
debris disks. The accretion of gas onto a giant planet requires a large solid core to form first, 
favoring a higher metallicity disk, whereas dust emission indicates the presence of smaller 
planetesimals that might be able to form in all disk environments. Another explanation may 
be that debris disks in high metallicity systems initiall y contain more material , but that over 



time all disks grind down toward similar masses (e.g. iDominik &: Decinl 120031 ) . In this case 
the detection of strong IR emiss ion is a reflection of a recent stochastic collision, rather than 
the disk's initial conditions (see 



Bryden et al.l (l2006al ) for further discussion). 



5.2. Age 



Collisions in a debris disk continually grind down the larger planetesimals, while the 
smallest dust can be removed by Poynting- Robertson drag and radiation pressure. One would 
assume that the overall disk mass must decline with time and, as expected, a correlation 
between stellar age and IR excess is observed, with debris disks more commonly identified 
around young er stars. While stu d ies concentrating on stars younger than 1 Gyr find a 



strong trend (jSpangler et al.l 12001 



the correlation is relatively weak (IBryden et al. 



Rieke et al. 



2005), among nearby solar-type field stars 



2006al ). In both cases, the evolution of the 



dust does not appear to be a steady decline. Observations of A stars find an overall decline 
in the average amount of 24 /zm excess emission on a ~150 Myr time scale, but the large 
variations on top of this trend suggest that sporadic collisional events a re able to dramati cally 
increase the amount of dust even at late stages in the disk's evolution (IRieke et al.ll2005j). As 



Decin et al. 


2000; 


Brvden et al. 


2006a) 



2001 



Figure [2] shows the resultant histogram of stellar ages. The ages for our main sequence 
stars are difficult to determine, with uncertainties in many cases of at least a factor of two. 
Where possible, we us e ages based on Ca II H&K line emission from the large compilation 
of IWright et al.l (120041 ) . Otherwise an average of values fou nd in the literature i s used. If the 
star was inferred to be young due to kinematic properties (IMontes et al.l 120011 ). we adopted 
that age. Table [1] lists the age data for each star. Although our target selection criteria do 
not explicitly discriminate based on stellar age, young stars (ages less than 1 Gyr) are not 
well represented in our sample due to their infrequent occurrence within ~25 pc of the Sun. 
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As in our earlier survey of nearby main- sequence stars (IBryden et al.ll2006al ). the stars 
with excess in this survey (marked with arrows in Figure [2]) have a weak but noticeable 
correlation with stellar age. No stars older than 7 Gyr have a significant amount of excess 
emission. The average age of stars with IR excess is 4.0±0.6 Gyr, compared to 5.6±0.4 Gyr 
for the sample as a whole. As discu ssed in the nex t sectio n, these trends are present in the 
combination of this sample with the IBryden et all (j2006al ) data. 



5.3. Spectral Type 

Observations of the general characteristics of debris disk as a function of spectral type 
are potentially a powerful tool for understanding the physical mechanisms responsible for the 
evolution of debris disks. The disk properties should be directly related to the stellar mass 
and luminosity in several ways. The mass of the protostellar disk from which the debris 
formed, for example, probably depends on the parent star's mass, as does its dynamical 
time scale. The stellar luminosity, though, is undoubtedly more important for debris disk 
characteristics, exerting a strong influence on the typical particle size (rbi OW out oc L*) as 
well as its temperature (T dnst oc L°" 2-0 ' 25 ). There are also observational biases linked to the 
brightness of the star, with cool dust seemingly easier to distinguish around hotter stars. 

The minimum L dnst /L+ based on the 70 /im flux, for example, is strongly dependent 
of stellar temperature (in eq. [3], detectable L dnst /L± is proportional to T~ 3 ). Hotter stars 
emit a lower fraction of their luminosity at infrared wavelengths, allowing for better contrast 
at those wavelengths. But while equation ([3]) is an observationally well-defined quantity, 
it contains no knowledge of the underlying disk physics. Naively, it appears to dictate a 
strong relationship between detectability and spectral type, i.e. it is easier to detect dust 
around hotter stars, but this may be misleading. For lack of any other information, the 
equation assumes that the dust emission peaks at 70 jum, thereby measuring the minimum 
L dust /L+. This assumed SED shape corresponds to a fixed dust temperature of ~50 K for all 
disks. One can instead consider disk models with a more physically motivated dep e ndenc e 



on spectral type. Instead of assuming a constant dust temperature, lHabing et all (120011 ). 
for example, assume the same dust location for all disks; in their models, the dust resides at 
50 AU independent of spectral type. In this case, dust temperature decreases with T*. In 
contrast with a simple reading of equation fl3]), the Habing et al. models have L dust /L+ more 
or less directly proportional to F dust /F+ for stars of type G and earlier. As in equation fl3j), 
lower stellar temperature makes it more difficult to detect dust emission relative to the stellar 
photosphere, but in Habing et al.'s models this difficulty is offset by a lower dust temperature 
for cooler stars, increasing the dust's 70 fim emission. 
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In Figure [U, the spectral types with IR excess stars are flagged with vertical arrows. 
A clear trend is readily apparent, with excess more frequently detected around earlier type 
stars. The detection rate drops from nearly 30% for the earliest type stars down to 0% 
for M stars. In fact, no stars with spectral type later than K0 are found to have excess 
emission (a sample of 23 stars without excess). This is consistent with previous survey 
result s that considered on ly part of the spectral range covered here. Our survey of F5-K5 
stars (IBryden et al.ll2006al ) found a detection rate of 13% within this limited spectral range, 
while a sample of ^ 30 images of nearby M stars yielded none with IR excess at 70 fim 
jGautier et al.ll2006h . 



A possible interpretation of the trend with spectral type is that it simply reflects the 
known correlation with stellar age. Earlier spectral type stars tend to be younger. Figure [TU1 
combine s information on spect ral type and age into a single plot for stars in this survey and 
those of IBryden et al.l (j2006al ). The trends previously identified are apparent - an upper 
limit to the ages of stars with excesses (filled symbols) of about 6 Gyr and a tendency for 
earlier type stars to have excess more frequently than later types. While the earliest type 
stars (F0-F3) are clearly younger on average, there is no clear evidence within the bulk of 
the sample that higher mass stars have more frequent excess because of they have younger 
ages. The formal correlation of excess with spectral type is even stronger than the correlation 
with age (correlation coefficients are -0.20±0.08. and -0.15±0.08 for spectral type and age 
respectively), further suggesting that spectral type is an independent indicator for IR excess. 
Unfortunately many of the latest type stars lack reliable age indicators, making it difficult 
to make any stronger conclusions. 



5.3.1. Comparison with Early Type Stars 



The detection rate of 70 jum excess around A stars is 33±4% (ISu et al.ll2006l ). more than 
twice that for the stars considered in this paper (13±3%). However, the A star and FGK 
star samples differ in both mass and age. We first consider the possibility that the different 
detection rates simply reflect an age evolution, rather than a spectral type dependence. For 
example, the youngest FGK stars have a detection rate somewhat higher than that within 
the sampl e as a whole: c onsider ing only systems with ages of 0.1-1 Gyr (and including the 
stars from IBryden et al.l (J2006a|)), 5 out of 19 young FGK stars have excess 70 /zm emission 
(= 26±12%). Similarly, the 70 fim excess frequency a mong the A sta rs drops with stellar 
age down to just 21±6% for A stars 0.3 - 1 Gy r old ( Su et all 2006 ). It is important to 
note, however, that many of the ISu et al.l (120061 ) observations are less sensitive than those 
presented here, relative to the stellar photosphere. Thus, their A star detection rate should 



-19- 



be regarded as a lower limit. Although the FGK and A star samples have stellar age as 
the most important correlating factor for IR excess, we cannot rule out a weaker but still 
important dependence of IR excess on some factor related to stellar mass such as luminosity 
or disk mass. 



5.3.2. Comparison with Late Type Stars 



Combining the observations presented here with those of iBryden et all (l2006al ) and 



Gautier et al.1 (120061 ). we have a total sample of 61 K1-M6 stars with no evidence of excess 
emission at 70 /im. Even considering only those stars whose photospheres are detected at 
70 /im with S/N>3 (42 of the 61 stars), this lack of excess detections is >3-a inconsistent 
with the ~15% detection rate around F and G type stars. As implied by Eqn (3), the 
contrast of dust relative to photosphere is, however, poorer for cooler stars which emit more 
of their energy in the infrared than hotter stars. The average upper limit to L^ ust /L^ for 
the 16 stars Kl or later with detected photospheres but no excesses in the SIM/TPF sample 
(S/N (70 fan) > 3 and X70 < 3) is Ld us t/£* < 9 x 10~ 6 compared with the average upper limit 
for 51 hotter stars with detected photospheres but no excesses, Ld U st/£* < 4 x 10~ 6 . Thus, 
one explanation for the lack of excesses around later-type stars is simply that the effective 
observational limits are a factor of two higher for the cooler stars. While observational 
selection effects make detection of IR excess around late type stars more difficult, the strength 
of this trend suggests that other explanations are needed. 

Another ambiguity in interpreting the correlation of excess with spectral type results 
from our limited knowledge of the location of the dust. If dust around later type stars is 
very distant from its central star, it will be too cool for detection at 70 /im. Figure [TT] shows 
how the dust temperature varies as a function o f spectral type for stars with excess from 



both this sample and from IBryden et al.1 ( 2006af ). We can only derive a dust temperature 



in the limited number of cases where we have a measured excess at both 24 and 70 /im 
(solid points). Otherwise, only upper limits can be obtained (Table Hj). For unresolved 
disk observations, the dust location cannot be determined without some knowledge of the 
underlying dust emission properties. Smaller dust with low emissivity can be just as hot 
as larger grains closer to the central star. Lines of constant orbital radius are shown in 
Figure [TT] under the assumption of either large blackbody grains (e = 1) or of small grains 
with emissivity = 0.01. Although the observed temperatures range from 80 to 170 K, they 
are all more or less consistent with emission from similar orbital locations: 10 AU for e = 1 
or 100 AU for e = 0.01. By implication, one would expect that dust around later K stars 
might have typical temperatures of ~50 K, ideal for detection at 70 /im, though none were 



-20- 



detected. 



Additional information on the location of the dust comes from IRS observations of F, 



G, and K stars with excesses (IBeichman et al.ll2006l ) which reveal that in almost all cases 
the inner boundary of the emitting region occurs at or interior to 10 AU. A theoretical basis 
for this preferred location of a radial distance of a few AU comes from the suggestion that 
the water-ice sublimation distance, or the "snowline" where the temperature falls below 170 
K, should m ark the onset of the region of giant planet formation and its remnants in the 
Kuiper Belt (jHayashilll98ll ; ISasselov fc Lecarll2000l ; iGaraud fc Lmll2006l ). Since the location 
of the snowline varies with stellar luminosity (oc L 25 ), there is no reason to expect a more 
distant, hidden reservoir of material unsampled by our observations around cooler stars. It 
is, of course, important to verify this expectation with obser yations at longer wav elengths 
such as MIPS 160 /im and in the sub-millimeter. Within the iGautier et al.l (120061 ) sample, 
for example, none of the 20 M stars examined at 160 /im show any excess emission, providing 
limits on Ldust/L* of 10 _5 -10~ 3 for material at ~50 AU. 

If the lack of debris disks around cool stars is real, then the dearth of material might 
reflect different formation mechanisms and evolutional history for the belts of planetesimals 
around low mass stars. Dust-producing collisions within these belts, for example, may 
require planetesimal stirri ng by larger, gas-gian t planets, whose frequency is thought to be 
lower for late- type stars (jLaughlin et al.l 12004 ). Alternately, the lack of IR excess might 
instead indicate a change in the physics of the smallest orbiting bodies as later type stars 
are considered, s uch as the increased relative importance of stellar winds in clearing dust 
from the system (IPlavchan et al.ll2005l ). 



6. Applicability to TPF 



The detection of other terr estrial planets is a long term goal for the astronomical com- 
munity (IMcKee fc Taylorll200ll ). NASA has spent considerable funds over the past decade 
on technology development and mission studies for a Terrestrial Planet Finder (TPF). One 
of the key astro-engineering issues revealed by those studies is the level of dust emission 
associated with target stars since exo-zodiac al emission is potentially an important source 
of photon shot noise (IBeichman et al.lll999l ). Thus, in addition to scientific interest, the 
incidence and distribution of material in the habitable zones, i.e. w here planets might h ave 
surface temperatures consistent with the presence of liquid water feasting et al.lll993l ). of 
nearby stars is of considerable technical importance. 
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6.1. Effect of Exo-Zodiacal Dust on Planet Finding 



As discussed in Appendix |AJ dust emission at the level of 10-20 times that of our 
own zodiacal cloud can impede planet searches (Figure [TBI) due to increased photon shot 
noise for either a coronagraph or an interferometer. Since this level is roughly 50-100 times 
less than that presently detectable with Spitzer, we can rule o ut only those stars with the 
most extreme zodiacal disks. Thus, HD 109085 and HD 69830 (IBeichman et al.ll2005al ) are 
unsuitable targets with strong excess shortward of 24 //m. However, the remaining stars in 
this sample and other samples pass the initial screening by having 24 /x m excesses, if any, les s 
than L dust /L* ~ 10~ 4 , corresponding to upper limits on \i EZ ~ 500 (IBryden et al.ll2006al ). 
Beyond these photometric constraints, IRS spectroscopy can push upper limits to factors of 
2-3 lower than MIPS alone and can also identify stars with small grain emission at 10 \xm 
( IBeichman et al.ll2006l ). 



In a few cases listed in Table El we can use the blackbodies fitted to the emission from 
the 5 stars with data at both 

24 and 70 fim (Table H|) to extrapolate the emission from this "Kuiper Belt" dust to 
the prime TPF-I wavelength of 10 /iin. The extrapolated emission is also given in units o f 
-^dust/At for material emitting at 10 //m ( cf. eq. [3] and eq. [21 of IBeichman et al.l (120061 )) 
relative to the solar system value of 10~ 7 ( iBackman fc Parescd Il993l ). Emission from any 
of this material located within the primary beam of the TPF-I telescopes (r <5 AU for 
a star at 10 pc observed with 3 m apertures) would be a noise source as described in the 
Appendix. However, this population of "cool" or "lukewarm" grains would not be located 
within a TPF-C pixel centered on the ~ 1 AU habitable zone and would not be a noise 
source at visible wavelengths. 

Unfortunately, however, the present observations cannot rule out an additional popula- 
tion of hotter grains located closer to the star which would either emit at 10 /xm or scatter 
in the visible. IRS observations in t he 8-14 jum region reach levels of just 1,000 times the 
zodiacal level (IBeichman et al.ll2006l ). It will take observations with nulling Interferometers 
such as the Keck and Large Binocular Telescope Interferometers which can spatially suppress 
the stellar component to measure directly the exo- zodiacal emission in the habitable zone at 
levels that could cause S/N or confusion problems for TPF. 

There is some cause for optimi sm, however. The "lu minosity function" of disks inferred 
from a variety of Spitzer samples (IBryden et al.l l2006al ). the rarity of extreme "hot" zodi- 
acal disks in the sam ple reported here and in other Spitzer papers (IBryden et al.l l2006al ; 
Beichman et al.l 120061 ) . and the apparent decline in the number of stars with excesses as a 
function of age (Fig. [2]) are all encouraging signs that the relatively clean example of our 
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solar system may be the norm rather than the exc eption. The ring-like st ructures seen in a 
number of resolved Spitzer disks, e.g. Fomalhaut ( Stapelfeldt et al. 12004) and e Eri (Back- 



man et al., in preparation) as well as in HST images (IKalas et al. 



200< 
200? 



20061 ) suggest that 



although the regions interior to the rings may not be completely empty due to a variety of 
mechanisms capable of transpo rting material inwa rd from the outer disk (comets, PR drag, 
interactions with planets, etc.; iHolmes et al.ll2002l ). these interior regions may have a quite 
low total optical depth, perhaps as low as the ~20% contribution inferred fo r material from 



Kuiper Belt material to the total amount seen at 1 AU in our solar system (ILandgraf et al. 
20021 ; bermott fc Kehoe! hooi iMoro-Martm fc MaJhot"rall2005l ). 



7. Summary 



We have searched for circumstellar dust around a sample of 88 F-M stars, by means of 
photometric measurements at 24 /im and 70 /im. We detected all the stars at 24 /im with 
high S/N and more than 80% of the stars at 70 /im with S/N > 3. Uncertainties in the Spitzer 
calibration and in the extrapolation of stellar photospheres to far-IR wavelengths limit our 
ability to detect IR excesses with 3-ct confidence to ~20% and ~50% of the photospheric 
levels at 24 and 70 /im, respectively. 

At these levels we have detected 12 of 88 ob jects with significant 70 /im excesses. Com- 



bined with an earlier study (IBryden et al.ll2006al ). we find an overall detection rate of 13 ±3% 
for mature cool stars. Beyond the single previously known 24 /mi excess within our sample, 
we detect two objects with 70 /im excesses and definite but weak 24 /im emission. Another 
two stars with 7 /im excesses h a ve 2 — a hints of 24 /im excesses. These results build 
on the finding of iBeichman et al.l (120061 ) that in many cases, objects with 70 /mi emission 
also had IRS spectra rising longward of 25 /tm to meet the 70 /im excess. These objects 
are all consistent with a disk architecture similar to our Kuiper Belt that is concentrated 
outside 5-10 AU. In this context we note that a number of the 70 /im sources are slightly, 
but significantly e xtended at 70 p. Th e detailed discussion of these objects is deferred to a 
subsequent paper (IBryden et al.ll2006bl ). The IR emission in these systems is different from 
the exceptional object H P 69830, which shows a disk architecture much more consistent with 
a massive asteroid belt (IBeichman et al.ll2005al ). 



Cross-correlating the detections of IR excess with stellar parameters we find no signifi- 
cant correlations in the incidence of excesses metallicity, but do find weak correlations with 
both stellar age and spectral type. The lack of correlation with metallicity contrasts with 
the known correlation between planet detections and stellar metallicity, and the expectation 
that higher metal content might result in a greater number of dust-producing planetesimals. 
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One significant finding is that the incidence of debris disks among mature stars is 
markedly lower for spectral ty pes later than KO tha n for earlier spectral ty pes. Combin- 
i ng da ta from this survey, the iBryden et al~ (j2006a ) F5-K5 survey, and the 



Gautier et al. 



( 120061 ) M star survey suggests an incidence of disks of 15±3% for F0-K0 stars and 0±4% 
for stars with types K2-M3. This lack of disks around later spectral types may be due to 
selection effects, lower initial disk mass, or different rates of dust creation or destruction. 

The disks that we are detecting have typical 70 /xm luminosities around 100 times that 
of the Kuiper Belt. If they also have inner asteroid belts 100 times brighter than our own, 
however, we would still not be able to detect this warm inner dust. The observed 70 jLtm 
excess systems could all be scaled-up replicas of the solar system's dust disk architecture, 
differing only in overall magnitude. These systems could have planets, asteroids and Kuiper 
Belt Objects as in our own system, but simply with a temporarily greater amount of dust due 
to a recent collisional event. Further observations of the warmer inner dust are necessary to 



addre ss this possibility Spitzer/IRS is particularly promising in this regard (IBeichman et al. 



20061 ) and is being pursued as part of a follow-up effort for some of the stars in this program. 



This publication makes use of data products from the Two-Micron All Sky Survey 
(2MASS), as well as from IPAC, SIMBAD, VIZIER, and the ROE Debris Disks Database 
website. We gratefully acknowledge the assistance of John Carpenter in reducing the IRAC 
data reported in this paper and we thank Angelle Tanner and Kate Su for helpful discussions. 
The Spitzer Space Telescope is operated by the Jet Propulsion Laboratory, California Insti- 
tute of Technology, under NASA contract 1407. Development of MIPS was funded by NASA 
through the Jet Propulsion Laboratory, subcontract 960785. Some of the research described 
in this publication was carried out at the Jet Propulsion Laboratory, California Institute of 
Technology, under a contract with the National Aeronautics and Space Administration. 



A. Noise Due to Exo-Zodiacal Emission 



In this section we make order of magnitude estimates of the impact of photon noise from 
exo- zodiacal emission on both visible light and mid-IR instruments (TPF-Coronagraph and 
TPF-Interferometer, respectively) designed to find neighboring planets. A detailed noise 
analysis of planet finding telescopes is b eyond the scope of this paper and the re ader is 
referred to other articles for more details (IBeichman fc Velusamylll999l ; lBrownll2005l ). 
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A.l. TPF-I, The Infrared Interferometer 

The use of a nulling inte rferometer to rej ect starlight and thereby reveal an orbiting 
planet dates to an article by iBracewelll (119781) and has been further investigated through 
studies of more sophisticated configurations flAngel &: Woolflll997l ; lLay et al.l 120051 ) . For a 
cryogenic system operating i n an orbit near 1 AU, the three dominant noise sources are 
( jBeichman fc Velusamylll999l . Table [5]): the stellar light that leaks past the interferometric 
null because of the finite diameter of the star, S^Leak] emission from the local zodiacal dust, 
Slz', and emission from the exo-zodiacal dust in the target star system that leaks past the 
interferometer, §EZ,Leak (see Figure [T41). At short wavelengths (<8 jttm), the stellar leak 
may dominate all other noise sources; longward of 20 /im emission from a 35 K telescope 
will become important; and at all wavelengths various systematic instrumental effects will 
be important. But over a broad range of wavelengths, the balance between S^Leak, Slz, 
Sez.Leofc controls the fundamental noise floor. Detector read noise and dark current can be 
ignored for broad band detection. 

In the background limit considered here, the total noise is given by the square root of 
the sum of all the individual photon fluxes reaching the detector. Rather than evaluate the 
absolute signal-to-noise ratio, S/N, we consider here the ratio of the S/N in the presence of 
exo-zodiacal emission, SNR(EZ), to the S/N in the absence of such emission, SNR(O): 
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In the above, S^Leak depends on the nulling configuration, the wavelength of operation 
and the angular size of the star. N ull depths of 10~ 5 to 10~ 6 have been demonstrated in 
the laboratory (Martin et al.l 120031 ) and for the purposes of this illustration, it suffices to 
take S^L eak = 10~ 5 F^. The emission from the local zodiacal cloud, Slz, is very complex 



in detail (IKelsall et al.lll998l ). but can be parameterized for our purposes as follows, Slz = 
t~lz B u (255 K)Q te i where B v is the Planck function, tlz is the vertical optical depth looking 
out from the mid-ecliptic plane at 1 AU, and Qtei is the diffraction limited solid angle of an 
individual telescope in the interferometer. A typical value of the zodiacal clo ud brightness 



towa rd the ecliptic pole from our mid-plane location is 12 MJy sr 1 at 12 /im (IKelsall et al. 
1998h . 



In the absence of more detailed information, the vertical optical depth of the exo-zodiacal 
dust in any system can be parameterized as a factor, Hez, times the Solar System's zodiacal 
dust. The emission from exo-zodiacal dust is then Sez{t) = 2^EZ T Lz( r )B u (T(r))Q te i, where 
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the factor of two accounts for the fact that in the exo- zodiacal case we are looking through the 
entire cloud and not from the vantage of the mid-plane as we do the local cloud. By analogy 
with the local zodiacal cloud (|Backmanl|l998l ) the vertical optical depth is assumed to fall off 
radially as r LZ {r) = T LZt iAu r Aii 3 - We also take T(r) = Tor AU as the equilibrium temperature 
for grains heated by stellar radiation and emitting in the infrared. Typical 1 AU values of 
(T n , ff) for large and small silicate grains are ( 255 K, -0.5) and (362 K, -0.4), respectively 
( iDraine fc Led 1 19841 ; iBackman fc Parescd Il993l ; iBeichman et al.l 120061 ) . The large and small 
grain brightness distributions are normalized to yield the same value at 1 AU. 

The effect of exo-zodiacal emission is modulated by the fringe pattern of the interfer- 
ometer which attenuates the bright central portion of the exo-zodiacal disk. To account for 
this effect we incorporate the fringe pattern of a particular nulling scheme ((9,<f)) where 9 
and 4> are the radial and azimuthal variables, respectively. In the simplified case of a face-on 
disk, the signal reaching the detector, $EZ,Leak is then given by the integral of Sez over the 
fringe pattern and the telescope solid angle: 
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for a star at a distanc e d. We adopt the fringe pa ttern ((9, </>) for the Dual Chopped Bracewell 
interfermoter (DCB; lLayll2004 lLay et al.ll2005l ) presently under study. Canceling out com- 
mon factors, the stellar leak term in equation flAlj) then becomes 
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To evaluate equation flA30 . we adopt a diffraction limited beam size of 
9 ma x = 0.6A/D = 0.5" for a D =3 m telescope at 12 fim. For a solar-type star at d =10 pc, 
the ratio of the exo-zodiacal contribution to that from the Solar System's own dust (eq. IA3I) 
is 0.06 hez ° r 0.24 \xez, for large and small grains respectively. Warmer, smaller grains fill 
more of the beam of the individual telescopes than the cooler, larger (blackbody) grains and 
thus contribute more noise. With this information in hand, Figure [TBI shows the variation of 
S/N as a function of exo-zodiacal brightness, hez, for two grain sizes. When the exo-zodiacal 
surface density hez is 10 times that of our solar system, corresponding to a 20-fold bright- 
ness increase, the S/N is reduced by a factor of ~2, necessitating an increase in integration 
time by a factor of ~4 to recover the original S/N. It is interesting to note the importance of 
grain size on this effect; the emission from the large grains is more centrally peaked and thus 
more effectively attenuated by the nulling interferometer than for the smaller grains which 
remain warm at quite large distances from the star. Since at least a few hours of integration 
time is needed to detect an Earth in the presence of a hez = 1 cloud, and days to carry out 
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a spectroscopic program ( Beichman 1998 ; Lay et al.l 2005 ). it is clear that studying systems 



with fi EZ > 10 — 20 will be difficult. 



A. 2. TPF-C, the Visible Light Coronagraph 

A similar analysis can be applied to an assessment of the effects of exo-zodiacal emission 
at visible wavelengths. There are some important differences however. First, the coronagraph 
takes in only the exo-zodiacal light from the immediate vicinity of the planet, not from the 
entire exo-zodiacal cloud (Fig. [TH right side). Second, the signal from an Earth (S p ), the 
residual starlight after a 10 -10 rejection ratio, and the local and exo-zodiacal signals are all 
more evenly balanced. Detector noise becomes a serious issue at medium spectral resolution 
(~75), but can be ignored in the broadband case. The analog of equation (lAlj) for the 
coronagraph becomes: 
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Since the local and exo-zodiacal emission enter the system through exactly the same 
solid angle, fl te i, the J^J term simplifies to 2^ EZ - For a planet 25 mag fainter than a V=4.5 
mag solar twi n at 10 pc, and assum ing a local zodiacal brightness of 0.1 MJy sr _1 at 0.55 



/im (Table El iBernstein et al.ll2002l ). we can evaluate the variation in S/N as a function of 
fiEZ- Figure [15] shows the decrease in S/N as the exo-zodiacal emission increases in the case 
of a face-on disk; an edge-on disk will increase the surface brightness and resultant noise. As 
with the interferometer, the effect of zodiacal emission in the target system is to lower the 
S/N by a factor of 2~3 at [Iez — 10. 

The relative effect of the exo-zodiacal emission is somewhat more pronounced for the 
TPF-C than for the TPF-I because the interferometer is dominated by the strong local zodi- 
acal background until very bright exo-zodiacal levels are observed. The intrinsic background 
level within the visible-light coronagraph is low (by assumption of an excellent 10~ 10 rejec- 
tion ratio) so that the exo-zodiacal emission more quickly plays a significant role in setting 
the system noise. 

A more detailed examination of the effects of the exo-zodiacal emission on the detectabil- 
ity of planets using TPF-C and TPF-I would yield absolute, not relative, sensitivity levels 
including the effects of disk inclination and confusion by structures, e.g. wakes and gaps, 
within the zodiacal cloud. These questions lie beyond the scope of this paper. 
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0.90 






1 


N 
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0.04 
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Average 
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HD 189245 
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5.65 
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HD 222237 


116745 


902 
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0.14 


E,M,N,T,V 
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HD 10360J 
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66 
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0.00 


E,T 


HD 38393 
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216A 
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7 Lep 
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0.30 


0.30 
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G3V 


5.56 


4.04 d 


12.90 






1 
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0.06 
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338A 






K7 


7.64 


3.99 
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73182 


570B 






M1V 


8.01 


3.90 e 




















a Age from Montes' paper, then from Wright's paper if available, otherwise an average of literature values 
b Known planet-bearing star 

c Star has a wide binary companion that was included in the survey, as listed at the bottom of the table 
d Star has one or more bad 2MASS values (err > 20%) 
c Star has JHK values from Johnson or other literature 



Note. — Spectral types from SIMBAD. Visual magnitudes are as quoted in SIMBAD, typically from the Hipparcos satellite; K magnitues are from 2MASS unless otherwise noted. 
References. — (B) iBarrvl l|l988l l: (B) iBarrvl l|l98Sh : (CS) ICavrel de Strobel et ail l|l996l . 120011 ): (C) IChen et al.l (|200lh ; (E) lEggenl l|l998h : (F) iFeltzing et aH l|200lt ): (I) llbukivama fc Arim 



Lambert fc Reddvl ll2004h; f La) lLachaume et"al] jl99S|); fM) ISheveley fc Marsakoyl Jl988l l; iMarsakov fc Sheyelevl (119951 s ) ; fMol iMontes et al.l [bOOll l; (N1 iNordstrom et al.l ll2004h ; (P) 
i"! ljl977h : (R'l lRandich etai\ l|l999h : fRPI iRocha-Pinto fc Maciell l|l998l l: (T1 lTavloiTl|2002l l: fWl lWright et al.l \200j ): (V) IValenti fc Fischej l|2005l l 



Table 2. Measured and predicted flux densities at 24 and 70 jum (in mJy) 



Star 


Fv,MIPS 


24 fim 

Fv,* 


F 'v.mips / 1 F V,* 


X24 


F V ,MIPS 


F v ,* 


70 ixva 
Fv.mipsI F v ,* 


S/N 


x b 70 


Ldustl L% 


GL 436 


38.7 


31.3 


1.24 


6.0 


4.0 ± 2.2 


3.5 


1.1 


1.9 


0.2 


< 5.4 x 10 -5 


GL 908 


96.3 


99.2 


0.97 


-0.7 


12.0 ± 3.3 


10.9 


1.1 


4.3 


0.3 


< 2.2 x 10" 5 


HD 739 
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142.5 


1.12 


1.4 


16.2 ± 3.8 


16.2 


1.0 


5.6 


0.0 


< 3.3 x 10~ 6 


HD 4391 


141.3 


150.7 


0.94 


-1.6 


19.5 ± 3.5 


17.2 


1.1 


8.4 


0.7 


< 3.9 x 10~ 6 


HD 4813 


191.8 


207.3 


0.93 


-0.9 


21.2 ± 4.6 


23.5 


0.9 


7.2 


-0.5 


< 2.7 x 10~ 6 


HD 10360 


247.7 


283.6 


0.87 


-1.6 


23.1 ± 6.1 


32.4 


0.7 


6.3 


-1.5 


< 4.6 x 10~ 6 


HD 16895 


492.7 


458.7 


1.07 


1.9 


51.2 ± 9.7 


51.9 


1.0 


8.8 


-0.1 


< 2.4 x 10~ 6 


HD 20794 


737.2 


770.8 


0.96 


-1.1 


94.3 ± 13.6 


87.9 


1.1 


28.7 


0.5 


< 1.2 x 10~ 6 


HD 22001 


233.7 


205.8 


1.14 


1.7 


25.2 ± 4.8 


23.3 


1.1 


7.8 


0.4 


< 2.7 x 10~ 6 


HD 23249 


2039.1 


1825.2 


1.12 


2.9 


207.0 ± 32.1 


206.5 


1.0 


24.5 


0.0 


< 1.5 x 10~ 6 


HD 23754 


383.3 


368.0 


1.04 


0.5 


46.6 ± 7.1 


41.8 


1.1 


13.8 


0.7 


< 1.5 x 10~ 6 


HD 25998 d ' e 


147.1 


128.5 


1.14 


3.6 


61.9 ± 5.7 


14.5 


4.3 


11.7 


8.3 


2.7 x 10~ 5 


HD 28343 


85.4 


98.4 


0.87 


-3.3 


-9.5 ± 14.8 


11.4 


-0.8 


-0.6 


-1.4 


< 8.8 x 10~ 5 


HD 32147 


224.3 


240.6 


0.93 


-0.8 


23.1 ± 6.3 


27.5 


0.8 


4.9 


-0.7 


< 7.6 x 10~ 6 


HD 36395 


250.7 


288.9 


0.87 


-3.3 


26.7 ± 17.6 


30.2 


0.9 


1.6 


-0.2 


< 4.9 x 10~ 5 


HD 38392 


203.7 


156.9 


1.30 


3.7 


18.8 ± 8.2 


17.7 


1.1 


2.4 


0.1 


< 1.8 x 10~ 5 


HD 38858 e -f 


131.3 


131.3 


1.00 


0.0 


153.7 ± 9.8 


15.0 


10.3 


16.0 


14.1 


1.0 x 10~ 4 


HD 39587 


483.6 


488.5 


0.99 


-0.3 


35.3 ± 12.1 


55.4 


0.6 


4.0 


-1.7 


< 4.0 x 10~ fi 


HD 40136 d - e 


553.3 


490.5 


1.13 


3.2 


90.7 ± 9.4 


55.6 


1.6 


21.3 


3.7 


4.0 x 10~ 6 


HD 46588 


150.8 


161.8 


0.93 


-1.7 


14.6 ± 3.7 


18.3 


0.8 


6.1 


-1.0 


< 3.0 x 10~ 6 


HD 48682 e -f 


188.4 


188.8 


1.00 


0.0 


256.9 ± 6.8 


21.3 


12.1 


46.3 


36.8 


1.1 x 10~ 4 


HD 50281 


175.8 


174.5 


1.01 


0.2 


17.1 ± 5.4 


19.9 


0.9 


3.8 


-0.5 


< 10.0 x 10~ 6 


HD 53706 


72.5 


73.4 


0.99 


-0.3 


7.0 ± 2.7 


8.3 


0.8 


3.0 


-0.5 


< 1.0 x 10~ 5 


HD 55892 


263.8 


243.0 


1.09 


2.1 


25.5 ± 5.7 


27.6 


0.9 


6.6 


-0.4 


< 2.2 x 10~ 6 


HD 62644 


370.6 


369.0 


1.00 


0.1 


46.2 ± 8.5 


41.9 


1.1 


8.2 


0.5 


< 4.0 x 10~ 6 


HD 63077 


226.4 


240.9 


0.94 


-1.5 


15.6 ± 8.0 


27.6 


0.6 


2.3 


-1.5 


< 6.2 x 10~ 6 


HD 67228 


215.6 


208.3 


1.04 


0.9 


14.4 ± 5.0 


23.5 


0.6 


4.0 


-1.8 


< 4.1 x 10~ 6 


HD 68146 


132.6 


132.5 


1.00 


0.0 


18.3 ± 3.6 


15.1 


1.2 


6.6 


0.9 


< 4.0 x 10~ 6 


HD 71243 


451.1 


434.5 


1.04 


0.5 


52.9 ± 8.2 


49.3 


1.1 


14.5 


0.4 


< 1.5 x 10~ 6 


HD 72673 


119.2 


132.8 


0.90 


-2.6 


9.2 ± 3.0 


15.2 


0.6 


4.6 


-2.0 


< 4.8 x 10~ 6 


HD 76653 


112.3 


109.1 


1.03 


0.7 


33.9 ± 10.5 


12.4 


2.7 


3.3 


2.1 


< 1.7 x 10~ 5 


HD 76932 


130.2 


140.7 


0.93 


-1.9 


15.2 ± 3.8 


16.1 


0.9 


5.3 


-0.2 


< 3.9 x 10~ 6 


HD 78366 


108.3 


115.3 


0.94 


-1.5 


16.8 ± 4.2 


13.1 


1.3 


4.6 


0.9 


< 6.8 x 10~ 6 


HD 79211 


195.4 


187.1 


1.04 


1.1 


18.0 ± 4.0 


21.4 


0.8 


7.4 


-0.8 


< 4.7 x 10~ 6 



Table 2 — Continued 



Star 


Fv.MIPS 


24 (im 


Fv.MIPS 1 F v ,* 


X24 


F„,MIPS 


F v> * 


70 

Fv,mipsI F v ,* 


S/N 


X70 


Ldust/ L% 


HD 81937 


563.9 


558.4 


1.01 


0.2 


70.0 ± 10.4 


63.3 


1.1 


16.6 


0.6 


< 1.0 x 10~ 6 


HD 81997 


287.8 


264.4 


1.09 


1.1 


28.2 ± 5.7 


30.0 


0.9 


8.0 


-0.3 


< 2.4 x 10~ 6 


HD 85512 


95.9 


110.9 


0.86 


-3.4 


11.0 ± 3.0 


12.8 


0.9 


4.7 


-0.6 


< 1.0 x 10~ 5 


HD 89449 


240.8 


255.1 


0.94 


-1.4 


24.6 ± 5.3 


28.9 


0.8 


8.2 


-0.8 


< 2.1 x 10~ 6 


HD 90089 e 


147.0 


146.7 


1.00 


0.0 


38.2 ± 3.6 


16.7 


2.3 


14.9 


6.0 


8.5 x 10~ fi 


HD 90589 


422.9 


409.8 


1.03 


0.8 


53.7 ± 8.5 


46.4 


1.2 


11.0 


0.9 


< 1.8 x 10~ 6 


HD 91324 


263.2 


273.9 


0.96 


-0.5 


50.5 ± 9.2 


31.2 


1.6 


6.4 


2.1 


< 5.3 x 10~ 6 


HD 100623 


183.8 


193.4 


0.95 


-1.2 


21.8 ± 4.1 


22.1 


1.0 


9.0 


-0.1 


< 4.0 x 10~ 6 


HD 102365 


353.9 


353.4 


1.00 


0.0 


34.0 ± 8.5 


40.2 


0.8 


5.7 


-0.7 


< 4.1 x 10~ 6 


HD 103095 


128.5 


133.8 


0.96 


-1.0 


9.3 ± 3.1 


15.4 


0.6 


4.4 


-1.9 


< 4.4 x 10~ 6 


HD 105211 e >-'' 


367.9 


363.4 


1.01 


0.2 


473.7 ± 19.8 


41.4 


11.4 


25.2 


21.8 


6.9 x 10~ 5 


HD 105452 


396.9 


409.1 


0.97 


-0.7 


43.3 ± 8.4 


46.6 


0.9 


9.3 


-0.4 


< 1.6 x 10~ 6 


HD 109085 d ' e - / 


589.2 


296.8 


1.99 


24.6 


198.2 ± 6.8 


33.5 


5.9 


42.7 


24.0 


3.3 x 10~ 5 


HD 129502 


536.0 


521.1 


1.03 


0.7 


51.8 ± 9.7 


59.1 


0.9 


13.1 


-0.8 


< 1.2 x 10~ 6 


HD 131977 


427.8 


397.4 


1.08 


1.9 


37.6 ± 7.4 


45.3 


0.8 


12.6 


-1.0 


< 3.2 x 10~ 6 


HD 132254 


127.5 


131.5 


0.97 


-0.8 


23.6 ± 3.4 


14.9 


1.6 


9.1 


2.5 


< 3.8 x 10~ 6 


HD 136352 


169.8 


175.1 


0.97 


-0.8 


17.5 ± 5.1 


19.9 


0.9 


4.2 


-0.5 


< 5.5 x 10~ 6 


HD 139664 e >-'' 


275.9 


251.9 


1.10 


1.2 


503.7 ± 9.2 


28.6 


17.6 


62.2 


51.9 


1.3 x 10~ 4 


HD 142267 


111.4 


107.5 


1.04 


0.9 


9.5 ± 2.7 


12.2 


0.8 


4.7 


-1.0 


< 4.2 x 10~ 6 


HD 147513 


201.1 


195.4 


1.03 


0.7 


17.3 ± 12.1 


22.1 


0.8 


1.5 


-0.4 


< 1.4 x 10~ 5 


HD 151288 


102.7 


115.9 


0.89 


-2.8 


13.1 ± 2.7 


13.4 


1.0 


7.4 


-0.1 


< 9.0 x 10~ 6 


HD 154363 


95.3 


110.5 


0.86 


-3.4 


9.7 ± 4.2 


12.7 


0.8 


2.6 


-0.7 


< 1.6 x 10~ 5 


HD 156026 


304.2 


306.9 


0.99 


-0.1 


30.6 ± 10.9 


35.2 


0.9 


3.2 


-0.4 


< 1.5 x 10~ 5 


HD 157881 


171.3 


192.8 


0.89 


-2.8 


13.1 ± 5.5 


22.2 


0.6 


3.0 


-1.7 


< 1.3 x 10~ 5 


HD 158633 e 


112.3 


125.1 


0.90 


-2.6 


56.7 ± 3.6 


14.4 


3.9 


19.9 


11.8 


4.1 x 10~ 5 


HD 160032 


243.8 


229.4 


1.06 


0.8 


44.3 ± 6.1 


26.1 


1.7 


9.4 


3.0 


< 3.3 x 10~ 6 


HD 164259 


246.4 


226.0 


1.09 


1.1 


8.3 ± 7.0 


25.6 


0.3 


1.4 


-2.5 


< 4.1 x 10~ 6 


HD 165499 


174.4 


182.1 


0.96 


-0.5 


15.9 ± 4.5 


20.6 


0.8 


4.9 


-1.1 


< 4.1 x 10~ 6 


HD 172051 


150.2 


147.4 


1.02 


0.5 


25.5 ± 11.5 


16.7 


1.5 


2.3 


0.8 


< 1.9 x 10~ 5 


HD 177565 


110.9 


109.8 


1.01 


0.2 


16.4 ± 4.9 


12.4 


1.3 


3.6 


0.8 


< 1.2 x 10~ 5 


HD 180617 


138.4 


111.6 


1.24 


6.0 


5.9 ± 13.4 


12.7 


0.5 


0.4 


-0.5 


< 9.6 x 10- 5 


HD 182488 


109.4 


108.7 


1.01 


0.2 


2.2 ± 10.7 


12.3 


0.2 


0.2 


-0.9 


< 2.7 x 10~ 5 


HD 185395 


284.9 


254.5 


1.12 


1.5 


34.0 ± 5.3 


28.8 


1.2 


11.1 


1.0 


< 2.0 x 10~ 6 


HD 187691 


212.0 


211.5 


1.00 


0.0 


21.0 ± 6.8 


23.9 


0.9 


3.6 


-0.4 


< 5.6 x 10~ 6 



Table 2 — Continued 







24 /^m 










70 /^m 








otcir 


Fv,MIPS 


f v.* 


r v,MIPSI r v,* 


X 2 4 


Fv.MIPS 


r v,* 


r v,MIPSI * v,* 


Q /1\T 


X70 


T I J c 


HD 189245 


120.8 


116.5 


1.04 


0.9 


6.6 ± 3.9 


13.3 


0.5 


1.9 


-1.7 


< 5.6 x 10~ 6 


HD 190406 


126.6 


133.0 


0.95 


-1.2 


22.6 ± 5.1 


15.1 


1.5 


5.0 


1.5 


< 7.9 x 10~ 6 


HD 191849 


179.0 


192.8 


0.93 


-1.8 


26.9 ± 4.1 


21.2 


1.3 


10.3 


1.4 


< 9.8 x 10~ 6 


HD 192310 


250.4 


236.4 


1.06 


0.7 


19.1 ± 5.6 


26.8 


0.7 


5.0 


-1.4 


< 6.3 x lO" 6 


HD 196877 


54.1 


59.4 


0.91 


-2.2 


2.2 ± 2.2 


6.8 


0.3 


1.1 


-2.1 


< 2.0 x 10" 5 


HD 198149 


2444.6 


2418.4 


1.01 


0.3 


254.6 ± 42.0 


276.7 


0.9 


38.5 


-0.5 


< 8.6 x 10~ 7 


HD 199260 d ' e 


120.1 


108.5 


1.11 


2.7 


42.8 ± 4.1 


12.3 


3.5 


11.7 


7.4 


2.1 x 10 -5 


HD 213845 


158.7 


144.7 


1.10 


1.2 


23.9 ± 4.0 


16.4 


1.5 


7.5 


1.9 


< 4.3 x 10~ 6 


HD 215648 


505.7 


516.7 


0.98 


-0.3 


47.4 ± 9.6 


58.8 


0.8 


12.7 


-1.2 


< 1.4 x 10~ 6 


HD 217357 


133.4 


139.3 


0.96 


-1.1 


18.6 ± 3.5 


16.1 


1.2 


7.5 


0.7 


< 8.6 x 10" 6 


HD 219482 d ' e 


140.9 


131.0 


1.08 


1.9 


65.4 ± 3.7 


14.9 


4.4 


22.5 


13.8 


2.8 x 10" 5 


HD 219623 e 


144.6 


139.7 


1.04 


0.9 


48.0 ± 3.8 


15.8 


3.0 


15.9 


8.4 


1.7 x 10~ 5 


HD 222237 


108.4 


112.0 


0.97 


-0.8 


20.5 ± 2.8 


12.8 


1.6 


9.9 


2.7 


< 7.1 x 10~ 6 


HD 265866 


84.8 


95.2 


0.89 


-2.7 




9.9 










Companions 




HD 10360J 


241.0 


337.7 


0.71 


-3.6 


23.8 ± 6.8 


38.2 


0.6 


6.5 


-2.1 


< 3.5 x 10- 6 


HD 38393 


765.4 


761.6 


1.00 


0.1 


57.3 ± 15.1 


86.4 


0.7 


7.3 


-1.9 


< 2.0 x 10~ e 


HD 50281B 


55.3 


43.7 


1.26 


6.6 


10.1 ± 4.6 


5.0 


2.0 


2.2 


1.1 


< 8.3 x 10~ 5 


HD 53705 


166.4 


192.1 


0.87 


-3.3 


17.3 ± 4.0 


21.9 


0.8 


7.3 


-1.1 


< 3.0 x 10- 6 


HD 79210 


191.4 


260.3 


0.74 


-6.6 


21.4 ± 5.1 


29.8 


0.7 


8.8 


-1.6 


< 5.6 x 10~ 6 


HD 131976 


274.4 


246.8 


1.11 


2.8 


15.6 ± 5.2 


28.2 


0.6 


5.2 


-2.4 


< 9.2 x 10~ 6 



a Significance of 24 fim excess (eq. pQ) 
b Significance of 70 excess (eq. [2]) 
c Minimum L^nat/L* from 70 )im emission (eq. [3]) 
d Star with excess 24 (im emission 
c Star with excess 70 ^m emission 
Star with resolved 70 /^m emission 



Table 3. IRAC Observations of M Stars" 





Spectral 


F ! ,(3.55 H 


F U (AA9 H 


F^(5.73 jum) 


F^(7.87 jum) 


Star 


Type 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


GL908 


Ml 


3.29 


2.22 


1.43 


0.82 


HD 36395 


M1.5V 


7.64 


5.62 


3.74 


2.22 


HD 191849 


K7/M0 


6.46 


4.12 


2.70 


1.58 


HD 265866 


M3.5 


2.75 


1.83 


1.22 


0.70 



a Flux density uncertainties are dominated by calibration uncertainties, typically 5% for 
the IRAC bands. 



Table 4. Dust Emission at 24 and 70 /im (in mJy) 



Star 


■fi/(24)d ua f 


F u (70) 


a 

dust 


rpb 
dust 


-^dust/ LI 


HD 25998 


18.6± 5.9 


51.1 


± 


9.6 


96±5 


4.5 


xlO" 


-5 


HD 38858 d 


< 15 (3a) 


193 


± 


25 


< 70 


12 


xlO" 


5 


HD 40136 


63± 22 


37.9 


± 


11 


165^ 


1.9 


xlO' 


-5 


HD 48682 d 


< 22 (3a) 


290 


± 


38 


< 68 


1.1 


xlO' 


-4 


HD 90089 


< 18 (3a) 


23.2 


± 


5.0 


< 120 


0.8 


xlO- 


5 


HD 105211 d 


< 45 (3a) 


521 


± 


73 


< 70 


7.3 


xlO" 


-5 


HD 109085 d 


292± 24 


212 


± 


28 


150±10 


15 


xlO" 


5 


HD 139664 d 


< 75 (3a) 


523 


± 


77 


< 78 


12 


xlO" 


5 


HD 158633 


< 13 (3a) 


45.8 


± 


7.8 


< 90 


3.5 


xlO" 


-5 


HD 199260 


11. 6± 5.6 


32.9 


± 


6.4 


94±5 


3.3 


xlO~ 


-5 


HD 219482 


9.8± 5.6 


54.5 


± 


9.0 


81±3 


3.6 


xlO' 


-5 


HD 219623 


< 17 (3a) 


34.8 


± 


6.5 


< 104 


1.6 


xlO" 


-5 



a Dust fluxes at 70 /im have been color corrected by 8%. 

b Blackbody temperature based on either the 24 to 70 /xm flux density ratio or the 70 /im 
flux density plus a 3a upper limit at 24 /im. 

c If only 70 /im data are available, L dust /L+ is from Eqn[3J If 24 /im and 70 /im data are 
available, L dnst /L^ obtained from Eqn|H 

d For the resolved sources, the emission is fit with an extended Gaussian profile, resulting 
in measured dust fluxes ~10% higher than from the standard aperture photometry. 



Table 5. Parameters for TPF S/N Calculations 



Parameter 


TPF-C 


TPF-I 


Wavelength 


0.55 ^m 


12 /iUL 


telescope 


3.5x8m 


rour, 3 m 






on 75 m baseline 


Beam Half Width 


39x17 mas 


500 mas 


Beam Area, Qtei 


5 x 1(T 14 


1.8 x lO" 11 


Local Zodiacal Emission (Ilz) 


0.1 MJy sr- 1 


12 MJy sr" 1 


Zodiacal Flux Density {Iiz^tei) 


5 nJy 


220 /iJy 


Stellar Magnitude 


V=4.5 mag (60 Jy) 


[12]=3.0 mag (1.7 Jy) 


Stellar Rejection 


io- 10 


10~ 5 


Stellar Leakage Signal 


6 nJy 


17 /xJy 


Planet Brightness 


6 nJy 


0.3 /iJy 



for a Solar twin at 10 pc. 



Table 6. Predicted Dust Emission at 10 /im 



Star 




Exo-Zodi a 




(mJy) 


(Solar System = 1) 


HD 25998 


0.041 


1.7 


HD 40136 


5.1 


42 


HD 109085 


23 


325 


HD 199260 


0.024 


1.1 


HD 219482 


0.004 


0.2 



a -krhist /£* at 10 Mm in u nits of 10 7 , c orresponding roughly t o that of the Solar System 
(IBackman fc Parescdll993l ). See eq. [2] of iBeichman et al.l (120061 ). 
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Fig. 1. — Spectral type distribution for stars in this SIM/TPF sample. The spectral types 
of stars found to have 70 /iin excess are highlighted within the histogram (slant-hash) and 
are individually flagged with arrows at the top of the plot. The length of each arrow is an 
indicator of the strength of 70 /zm excess relative to the stellar photosphere. We find that 
70 /im excess is more readily detected around early type stars. 
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Fig. 2. — Age distribution for stars in this SIM/TPF sample. The ages of stars with 70 jum 
excess are highlighted within the histogram (slant-hash) and are individually flagged with 
arrows at the top of the plot. The length of each arrow is an indicator of the strength of 
70 /j,m excess. There is a weak correlation between the detection of IR excess and the stellar 
age, with no stars older than 7 Gyr having excess emission. 




Fig. 3. — Metallicity distribution for stars in this SIM/TPF sample. The ages of stars with 
70 (Ltm excess are highlighted within the histogram (slant-hash) and are individually flagged 
with arrows at the top of the plot. The length of each arrow is an indicator of the strength of 
70 (Ltm excess relative to the stellar photosphere. There is no correlation between metallicity 
and the detection of IR excess. 
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T 1 1 1 1 1 1 1 1 

Stars with 70/xm excess 




at 24/xm 

Fig. 4. — Distribution of 24 /iin fluxes relative to the expected photospheric values. A 
Gaussian distribution with 10% dispersion (solid curve) is shown for comparison. One star 
(HD 109085, a star previously identified as having excess emission) clearly stands out from 
the main population. The broad dispersion within this population is due to a variety of 
factors. Some stars have poor estimates of the stellar flux at 24 /zm due to poor near-IR 
data or photospheric models, particularly for the set of late K and M stars marked 'poor 
extrapolation' in the figure. The spread of values is also increased by sources with true, 
weak excesses (at the level of ~10% above the stellar photosphere). Stars with excesses at 
a longer wavelength (70 /im) are shown with black shading. 
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• This survey 

x Bryden et al. 2006 

a Gautier et al. 2006 



* 



i 



F0-F4 F5-F9 G0-G4 G5-G9 K0-K4 K5-K9 M0-M4 

Spectral Type 



Fig. 5. — Average 24 fim color relative to 2MASS K s band (2.16 /xm) as a function of spectral 
type. Stars with excess emission or with poor K s measurements are excluded. Error bars 
indicate the err or on the mean value within each bin (not th e overall dispersion) . Stellar 
colors from the iBryden et al.l (l2006al ) F5-K5 survey and the iGautier et al.l (120061 ) M star 
survey are also shown for comparison. The trend is relatively flat over most of the range, 
with significantly red colors only seen among the M type stars. 



-47- 



o 
M" 



O 

V) 

o 

■+-> 

o o 

CM 



~i 1 1 r 



~i 1 r 



1 2 3 4 5 



Excess (>3— a) 

No significant excess 



10 



15 



Fig. 6. — Distribution of 70 /im fluxes relative to the expected photospheric values. While 
most stars cluster around unity (consistent with emission from the star alone) many show a 
high degree of excess emission attributable to circumstellar dust. 



Fig. 7. — Composite image of the field surrounding HD 105211 (marked by a plus sign). In 
addition to our MIPS 24 /im (green) and 70 /im (red) images, the 2MASS K s band image is 
overlaid in blue. While dim background stars show up as blue points, the cool Mira variable 
star CL Cru (triangle) has strong 24 /xm emission with a much broader PSF, resulting in an 
overall green color. This bright star comtaminated the broad IRAS scanning beam (dashed 
yellow rectangle) prohibiting the detection of dust around HD 105211. In the MIPS image, 
HD 105211 is well resolved, showing strong excess emission at 70 /im. The neighoring star 
CPD-63 2145B (asterisk), while detected at 24 /im, does not give off significant emission at 
70 /im. 



-48- 




1 10 100 

A (//in) 



Fig. 8. — Spectral energy distributions for two stars with IR excess at both 24 and 70 /zm. 
The emission from the two stars has been scaled such that their photospheres overlap, empha- 
sizing the difference in far-IR emission. The observed fluxes at each wavelength are shown as 
open circles that are fit with a combination of emission from the stellar photosphere (dotted 
line) and from orbiting dust. The dust emission of HD 40316 is fit with 170 K dust, whereas 
HD 219482, with stronger 70 /zm and weaker 24 /zm emission, is fit with cooler, 80 K dust. 
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Fig. 9. — Area and dust mass estimates for stars with 70 fim excess em ission. In addition to 
our 12 stars with excess, the 7 excess stars from iBryden et al.l (l2006al ) are also shown. Dust 
temperatures (Table H|) are translated into orbital radii assuming either large blackbody 
grains (bottom axis) or small grains with emissivity 0.01 (top axis). Stars with excess 
measured at both 24 and 70 ixm are shown as solid points, while those with only upper 
limits for the dust temperature are shown as open circles. Error bars are added to each 
point based on the 1-er uncertainties in the dust temperature; for systems with upper limits 
on the dust temperature, an arrow is plotted with length/direction based on an assumed 
10% uncertainty in the temperature. Dust masses (right axis) are calculated assuming a 
typical grain size of 10 fim. Both dust area and mass are calculated under the assumption 
of blackbody grains (unity emissivity); for an emissivity of 0.01, both area and mass are 
a factor of 100 larger. The detection limits, which depend on the stellar temperature, are 
shown for a F0 star (dotted line) and for a K0 star (dashed line). 



-SO- 



LD 



-1 1 1 r - 



~i 1 r 



O 



A O 



o 



o 



-1 1 1 r - 



T 



n r~ 



o This survey 
a Bryden 2006 
filled: IR Excess 



o 



A 



O 



o o 




• o 




o • 8 




• • 


o 


o 






o 


1 


1 • 


F0 


F5 



o 

o 

• A 

• o 

A 

o A 

A 

$A 

A 



A 



A A 

A ° A 



A 



o 

A 



O 



A 

o 



A « 

£ A 



O 



▲ 

A 
_1_ 



A 



O 



A 
A 



o o 



A 
A 



o 

A 



A 



o 



. ° . 



GO G5 K0 

Spectral Type 



K5 



MO 



Fig. 10. — Stellar age as a function of spectral type for stars with known age s. The stars 
from this survey are marked with circles, while those from iBry den et al.l (j2Q06al ) are marked 
as triangles. In both cases, stars with IR excess are marked as filled symbols. 
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Fig. 11 . — Dust temperature for stars with IR excess in this sample and from lBryden et al.l 
(j2006al ). For stars with excess measured at both 24 and 70 /im (solid points), the dust SED 
is fit with a representative temperature. Those stars with a single measurement of excess at 
70 /im only have (3-cr) upper limits. With some assumptions for the grain properties, the 
dust temperatures can be translated to orbital distances. Several lines of constant distance 
are shown for comparison. The observed systems with measured dust temperatures are 
mostly consistent with large blackbodies orbiting at ~10 AU or small, low-emissivity grains 
at ~100 AU. There is no clear evidence for orbital distance changing as a function of spectral 
type. 
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Fig. 12. — Detection frequency of planets and IR exc ess as a function of ste llar metallicity. 
For the stars presented in this paper, plus those of iBryden et al.l (j2006al ). the detection 
rate of 70 /im excess emission is shown as open circles. The distribution has no trend in 
metallicity; all points are consistent with the average detection rate for the entire sample 
(14%; dotted line). This is in contrast to the dependence of the planet detection rate on 
metallicity for a similar sample of nearby stars (x marks), which Fischer &: Valentil (120051 ) 
fit with a metallicity-squared relationship (dotted line). 
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Fig. 13. — Distribution of r correlation coefficients for two series of monte carlo simulations. 
Stars were selected either completely randomly (left histogram) or proportionate to their 
metallicity squared (right histogram), the relationship observed for planet-bearing stars. 
The arrow at the top shows the value of r observed within our data, strongly inconsistent 
with the planet-metallicity relationship. 
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5 AU at 10 pc 



Visible Coronagraph 
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0.4x0.17 AU at 10 pc 



Fig. 14. — The interferometer (left) takes in all emission from a face-on exo-zodiacal disk 
that fits within the diffraction-limited beam of a single telescope. The intensity is, however, 
attenuated by the interferometric fringe pattern shown as vertical black bars. The corona- 
graph (right) takes in exo-zodiacal light only within the area of a single diffraction-limited 
pixel; the values shown here are appropriate for a 3.5x8 m telescope under consideration for 
TPF-C. 
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Fig. 15. — The effect of exo-zodiacal dust emission on TPF signal-to-noise ratio. The 
horizontal axis gives the vertical optical depth of the exo-zodiacal disk normalized to that 
of the Solar System (/jllz)- Note that a value of fi E z = 1 in a target system corresponds to 
twice the emission we see from our zodiacal cloud, e.g. using COBE or IRAS, since we view 
our cloud from its midplane. The upper two curves show the falloff in relative S/N for the 
interferometer as the amount of exo-zodiacal emission increases, with large and small grain 
sizes considered separately (solid and dotted lines respectively). The lower curve (dashed 
line) shows a similar trend for a coronagraph viewing a face-on disk. In each case, the 
signal-to-noise ratio is shown relative to observations of a system with no dust emission. 



